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A general scheme for devising efficient cluster dynamics proposed in a previous letter 
[Phys. Rev. Lett. 72, 1541 (1994)] is extensively discussed. In particular the strong connection 
among equilibrium properties of clusters and dynamic properties as the correlation time for magne- 
tization is emphasized. The general scheme is applied to a number of frustrated spin model and the 
results discussed. 



I. INTRODUCTION 



' The cluster formalism introduced by Kasteleyn and Fortuin (KF) Q and later developed by Coniglio and Klein 



(CK) for the ferromagnetic Ising model has greatly enhanced the understanding of critical phenomena in terms of 
geometrical concept. Moreover, based on such a formalism Swendsen and Wang (SW) Q have introduced a cluster 
dynamics which drastically reduces the critical slowing down in Monte Carlo (MC) simulation of ferromagnetic spin 
models. In the SW dynamics spins belonging to the same cluster are flipped in one step as opposed to the single 
spin dynamics, where spins are flipped one at time. The efficiency of the SW dynamics stems from the fact that in 
the KF-CK formalism the clusters represent correlated spins; therefore if one spin in a cluster is flipped all the other 
, spins in that cluster will successively tend to flip coherently. Consequently by flipping in one move all the spins in 
' the same cluster results in a much faster dynamics. The SW algorithm has been extended and applied efficiently to 
' several unfrustrated spin models. @-|^ 

Unfortunately, the SW dynamics based on the direct extension of the KF-CK cluster formalism to frustrated spin 
-p^ systems does not show any reduction of the relaxation times. The reason for that is due to the fact that now 

d the KF-CK clusters do not represent any longer correlated spins. Recently Kandel, Ben-Av and Domany have 

^ introduced a new type of MC cluster algorithm for the particular case of the fully frustrated Ising model which 
is able to reduce the critical slowing down. Attempts to use their algorithm to other frustrated models has been 
^ , satisfactory in few cases and discouraging in others A different approach based on the definition of 

^ quasi-frozen clusters in spin glasses looks very promising, but the implementation of a related cluster dynamics has 
't^ not been yet explored . More recently, based on the approach of Kandel et al., we have proposed a general criterion 
and a systematic procedure to define clusters and related efficient cluster dynamics for frustrated spin models | |20[ | . In 
particular, we have applied our general criterion to a class of fully frustrated Ising spin models on square lattices where 
5^ \ the relative strength between the interactions can be varied. For any value of the relative strength, without invoking 
"ad hoc" algorithms for each case, the same general criterion generate a Monte Carlo dynamics which dramatically 
reduce the critical slowing down. 

The aim of this paper is to illustrate in more details the criterion proposed in Rcf. [ pO[ and apply our method to a 
number of frustrated spin models. In section II, III we discuss the extension of cluster formalism to frustrated spin 
models. We stress that for the unfrustrated spin model the clusters percolate at a temperature Tp that coincides 
with the critical temperature Tc, while for the frustrated models Tp is larger than Tc- In section IV following the 
approach of Kandel et al. we introduce a large variety of cluster definitions which contains as a particular case the 
KF-CK clusters. We then illustrate a procedure to approach systematically, in successive order of approximations, 
a cluster definition for which the percolation temperature Tp becomes closer and closer to the critical temperature 
Tc- In section V, VI and VII we check our procedure by comparing percolation quantities pl|j2^ ] and thermodynamic 
quantities on a variety of frustrated models using MC simulations. For a number of frustrated models without 
disorder we find to second order that the clusters percolate at a temperature Tp numerically indistinguishable from 
the critical temperature Tc- For disordered frustrated models like Spin Glass (SG) the convergence of Tp towards the 
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SG critical temperature Tsg is very slow. In section VIII we implement a cluster MC dynamics based on the novel 
cluster definition. We show that the dynamics is very efficient with drastic reduction of the relaxation time for those 
frustrated systems, introduced in section V-VII, for which the cluster definition leads to Tp ~ Tc- In the appendix we 
show that the cluster dynamics fulfills detailed balance and briefly discuss the ergodicity problem. 



II. THE CLUSTER APPROACH IN FRUSTRATED SYSTEMS. 



It is well known that he partition function of a ferromagnetic Ising model can be written in terms of the clusters 
of an equivalent percolation model Rj2| . A similar result can be also obtained for Ising systems where frustration is 



present |12|. The aim of this section is to recall and discuss those results which will be useful in the following. 



Let us consider the Ising Hamiltonian 

mSr}) = - E Ji^^3S^S, - 1) (1) 
<ij> 

where = ±1 is the sign of the quenched interaction and J > the interaction modulus. The interaction configura- 
tion {tij} can be a deterministic periodic structure or a disordered one. Hamiltonian ([^) is said to contain frustration 
if at least one closed path C exists such that Y\(^^ ^^^^ = —1. 

Within the CK approach |^,^^ we introduce "bonds" between nearest neighbor (nn) spins satisfying the in- 
teraction probability p = 1 — e^'^^'^ . The weight for each given configuration of spin {Si} and bond C is given 
by 

WcK{{S,},C)^p\''\l-p)\^\ n ^s.s, n (l-'^s.s,) (2) 

<i,j>eCF <i,j>(£CA 

where Cp (Ca) is the subset of C which covers ferromagnetic (antiferromagnetic) bonds, i.e. bonds with = 1 
(e,y = —1), and p = 1 — e"'^^'' . The product takes into account the fact that bonds can link only spins which satisfy 
the interaction. The clusters are defined as maximal sets of spins connected by bonds. It can be shown [ p^ that in 
this case the partition function becomes 

^ = ^*p|c|(i_p)|A|2iV(c) (3) 
c 

where means that the sum is over the bond configurations which do not contain frustration. Furthermore, the 
presence of positive and negative interaction implies that the following relations can be proved 

(5.) = {1°^)CK ~ {iri)cK (4) 

and 

{S^Sj) = (7II )cK - {ll)cK (5) 

where (7if ) is 1 if the spin i is "up" ("down") and belongs to an infinite cluster, otherwise it is 0, and (7]^) 
is 1 if the spins i and j belong to the same cluster and are parallel (antiparallel), otherwise it is 0, and where (...) is 
the usual thermodynamic average for a fixed configuration of interaction {Jcij} and {■■■)ck is the average over spin 
and bond configurations weighted with (^). The percolation quantities are instead given by 

^ {^ncK^ (7n + {^^)CK (6) 
Pi] = (7u)cx= {ll)cK + {ll)cK (7) 

where 7°° = 7^^ + and 7^ — 7-^ + "ff- , and Pi is the probability that the spin at site i belongs to the 00 cluster 
and Pij is the probability that the spins at site i and j are in the same cluster. 
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It is clear that without frustration the relations valid in the ferromagnetic case are recovered with only trivial 
differences; in fact, in this case, eqs. (^, ||) become 

\{S^)\^{lr)cK (8) 

and 

\{S,S,)\^{j,,)cK. (9) 

Examples of models without frustration for which eqs. (^) and (||) are satisfied are the antiferromagnetic Ising 
model on a square lattice and an Ising model with interaction Jij — JuiGj where ct^ = ±1 are quenched variables. 
It is easy to realize that for any fixed configuration of Ui, although the interactions may be positive and negative, all 
the loops are unfrustrated. 

From (H) and it follows that unlikely the ferromagnetic case the critical temperature does not coincide with 
the percolation temperature. In fact, defining the critical temperature Tc as the temperature at which the Edwards- 
Anderson order parameter qEA = jj^i^i)"^ vanishes and the percolation temperature Tp as the temperature at 
which the percolation probability poo = jf^iPi vanishes, from eq. (|^) follows that Tc < Tp |l2|. (In the definition 
of qEA the bar stands for the average over the configurations of interaction {Jeij} and N is the number of spins. 
psf ) This result has been verified numerically for a number of different systems as 2d and 3d spin glass [^,|lO| , fully 
frustrated model |]lO| and frustrated XY model ||ll[]. 

Similarly from eq. correlation and connectivity do not coincide any more. In fact two spins instead of being in 
the same cluster may be parallel in one configuration and antiparallel in another (see Fig. 1). Although these two 
configurations will both contribute to the connectivity they will interfere and strongly reduce correlations. Therefore 
for T > Tc defining gij = (SiSj) it follows from (||) and (^) \gij \ < Pij. 



III. GENERALIZATION OF THE CLUSTER APPROACH. 

The main result of the previous section is that when frustration is present the KF-CK clusters do not represent 
anymore correlated spins. As a consequence the percolation temperature Tp is higher than the thermodynamic critical 
temperature T^. We will generalize now the KF approach ^ in order to define new clusters which represent correlated 
spins even in the frustrated case. The new clusters will reduce to the usual KF-CK clusters in the unfrustrated case. 

We will achieve this goal in two steps. First, in this section, following the approach introduced by Ben Av at 
al. fl^ , ^ (BKD), we will consider a large class of clusters which contain as a particular case the KF-CK clusters. 
Second, in the following section, wc will give a criterion to choose systematically the "right" clusters in successive 
order approximations in such a way that pair connectness and pair correlation functions tend to coincide. We will 
consider an approximation good enough when the percolation temperature Tp has approached the critical temperature 
Tc- 

Let us consider a square lattice with ferromagnetic and antiferromagnetic interactions |^7| and let us focus our 
attention on a single isolated plaquette. This can be either unfrustrated or frustrated. Following Ref. [ iQ], it 
is possible to generalize the KF procedure by assigning globally to each bond configurations on a given plaquette or 
even on larger spin blocks taken as "elementary units" (see Fig. 2)] its own bond probability. Within each elementary 
unit or plaquette each bond configuration probability is independent on each other, and thus, it is not given, as in 
the CK formalism, by the product of bond probabilities of single pairs of spins. 

We consider, to begin with a specific case, a checkerboard partition of the square lattice (see Fig. 3) and we take one 
of the two sets of plaquettes (the shaded or the unshaded ones) as the set of "elementary units" on which we make our 
independent choices. We proceed further by generalizing the KF approach. Consequently we "dilute" the couplings 
on the plaquette, replacing them with a set of new interaction configurations which contain only J' ^ oo or J' ~ 
interactions, We consider the generic hamiltonian (1^). In sec. V-VII we will consider specific examples. The 
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possible interaction configurations are shown in Fig. 4. We assign a bond weight wq to the interaction configuration 
Co of the first row (no bonds), we assign the same weight wi = wi...— to all the configurations of bonds, C1...C4 
in the second row, and so on. The symmetry of the plaquette allows us to choose the same weight for symmetric 
configurations, i.e. members of the same row in Fig. 4. 

The requirement of the equality between the partition function of the original model and that of the diluted model 
gives 

M 

<i,j>Gplaq. a—0 <i,j>Gc^ 

where the sum is over all the possible interaction configurations on the plaquette {M — 15 for the example in the 
Fig. 4). With < i,j >G Cq, we mean the n.n. spins connected by J' interactions in Ca and with < i,j >G plaq. the 
n.n. spins on the plaquette. The previous relation can be rewritten as 

M 
a=0 

where 

'Hplaq.m},Ca,) = - J2 J'i^vS^SJ-l) (12) 
<i,j>£Ca 

and 

npia,.i{S,}) = ~ J2 Ji^vS^S,-!) (13) 

<i,j>£plaq. 

are the energy of the plaquette viewed as "elementary unit" in the new interaction configuration and in the original 
system. 

Of course such procedure can be repeated for every Hamiltonian that can be written as sum of elementary unit 
energies, i.e. 

n{{s.})^J2^ias,}), (14) 

I 

whose block I we dilute and obtain the stochastic Hamiltonian 

7i({5,},C)-^7i,({5,},c„(z)) (15) 
/ 

with C = U, Ca(l). 

The equivalence between the original model and the diluted one is obtained imposing (cfr. eq.(|ll|) ), 

E ^^awe-^'^'^^^-^'^"''') = e-'^^^"^"" (16) 

where the Wa{i) is weight with which the configuration Ca(i) on the Zth elementary unit occurs. 
Furthermore, the partition function can be written as 

z=Y^Y[e-^^'^''y= y: n^"''^'"^'^""'"^^"^) ■ (17) 

{Si} ; {s,},c I 

Spins that are connected by infinite strength interaction are frozen while the others do not interact. Thus the 
dilution of the original Hamiltonian is also called a "freezing and deleting operation" . Of course we have: 
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,-,miiS.M = 5{5.},c„,, (18) 

where S^Si},Ca, 1 o^' depending whether or not the spm configurations satisfy aU the oo strength interactions in 
the interaction configuration Ca{l) of the Z-th plaquette. Two spins connected by an infinite strength interaction wiU 
be frozen in the configuration which satisfy the interaction. On the entire lattice we can define: 

hs,},c^Y[S[s.},c^^„= n ^s.s, n i^-^s.s,) (19) 
I <i,j>eCF <ij>eCA 

where Cp and Ca are defined hke in eq.(||). 
Let us observe that 

M 

n ^c(/)'5{sa,c = n K)"°^''^'5|5.},c (20) 

where na{C) is the number of elementary units on which we have chosen the ath configuration (with a ~ 0, ...M) in 
the given interaction configuration C (with na{C) — Mu where Mu is the number of elementary units). 
Therefore, from eqs. ([l7|), ( |l8| ) and (|20|), we obtain 

M M 

^ = E n(-")"°'''^ E hs.},c = E* n(-")""'''^2^^''^ (21) 

C a=0 {S.} C Q=0 

where stands for the sum over all the interaction configurations that are compatible with at the least one of the 
possible spin configurations (i.e. all the unfrustrated graphs). Eq. (21) is the generalization of eq.(||) which can be 
recovered considering a pair of n.n spins as elementary unit. 

Equation (^) can be also obtained following the CK approach where the clusters are defined in the original system 
introducing fictitious bonds between spins satisfying the interaction. Given a spin configuration {Si}, the probability 
to realize a configuration of bonds Cq,(j) on each unit / is given by 



where Wa satisfy eq.(^^. Due to (|16|) and (|l8|)these probabilities are normalized for any spin configuration 
li*^*}) ^ "^^^ Kronecker delta assures that the bonds are thrown only between spins satisfying 
the interaction. For the entire system the weight for a given configuration of spins {Si} and bond configuration 
C = U; Ca(0 is given by 

W^({50,C)=nnc„(,)|{5.})e-^«»^>»=n^"(0'5{S.},c„<o (23) 
I I 



where (14h have been taken into account. Finally from (h9h and (EQ) we have 



M 



W{{S.}, C) = n («'-)"" ''^^^{so.c. (24) 



Q=0 



Summing over the spin and bond configurations we recover eq. (pi|) . The advantage of this approach is to make 
clear that both spins and bonds can be defined in the original system where the clusters are defined as maximal 
sets of spins connected by bonds. To calculate the statistics of the CK clusters we have to generate equilibrium spin 
configurations first. Then, for each equilibrium configuration we can assign a bond configuration on each unit I with 
probability given by (p^. 

Summarizing, following the approach of BKD, we have defined a vast class of generalized percolation models 
equivalent to our original spin model eq. (p]) . When one reduces the elementary unit to a single spin pair one recovers 
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KF-CK solution and that for larger units it is always possible to find a solution in the form of product of KF probability 
p = 1 — e"'^^"' , or, generally, solutions that are factorization of probabilities for a subpartition of the elementary units. 

The generalization discussed does not solve automatically our problem of recovering the identity between cluster 
connectivity and spin correlation function, eq.(^), when frustration is present. However, the great freedom given by 
eq. (|l|) stiU gives hope to find solutions for w^f^i^, in such a way to achieve the equality @ at least in an approximate 



way. Then it is crucial to find a criterion to select among the many different possibilities offered by eq.(16), those for 



which (H) holds. This is the second step needed to achieve our aim and will be discussed in the next section. 

IV. CONDITIONS BETWEEN CORRELATION AND CONNECTIVITY. 



For each percolation model defined in the previous section (satisfying eq.(|l6|)), it is straightforward to generalize 
the relations (jj) and (||) provided that now the average over the spin and bond configurations has to be computed 
with weights given in eq. (Q) . In particular relation is also valid when we consider a sub-system made of a single 
unit I, namely 



{S^S,)i = (7II )/ - (7,^),, (25) 



for each i and j on the unit I. Here 



= E{5.|e-/3«^((s.}) (26) 

where the sum is over all possible spin and cluster configurations on the unit / and P{ca\{Si}) is given by (|2^). When 
the quantity to average is function only of the spin variables like {siSj)i, due to eq. ([22[), eq.(p6[) simplifies 

" E{s,}e-/3«'«s.}) (2^) 

Our aim is to find among the large class of solutions of eq.([l^) a solution for Wa{i) in such a way that the equality 
(^ is satisfied. In this way the percolation temperature Tp would coincide with Tc and the clusters so identified will 
be characterized by percolation critical exponents equal to thermodynamical critical exponents. Since this is not a 
trivial task, we seek solutions which fulfill approximately eq. (|^) by imposing the condition on a sub-system made of 
a single unit, namely 

\{S,S,)i\^{j,,)i. (28) 

To find the solution of (|l6|) and (|2^) may still be complicated due to the large number of unknown (2" where n is 
the number of edges in the unit cell). Technically the problem can be simplified if we make use of eq. (|2^) which is 
always valid if eqs ([l6| ) are satisfied. 

Eq. ( p8| ) together with eq. (p5| ) implies for each pair ij in the elementary unit 

{jbi - if {s.s,)i > 

(4)1 - if {S,S,h < (29) 
(7,j), - if {S^Sj)i = 0. 

If ( p9| ) and (^6|) are satisfied, from ( |2^ ) follows that also ( p^ ) is satisfied. As we will show in the next section in a 
specific example it is easier to impose (|l6| ) and (|2^) than to impose the equivalent conditions (^6|) and ( p8| ) . 

It is clear that the larger the unit the better eq.(||) is satisfied. It also clear that there is no guarantee on how 
good the different approximations are and how fast eq.(|p) is approached by increasing the unit size. However, this 
procedure, by increasing the size of the unit, allows a systematic way to improve the approximation. 
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Eqs.(^6|) and ( p8| ) introduce a set of independent conditions whose number depends on the size and the symmetry 



of the chosen elementary unit. In general we are not able to know if the conditions introduced by (16) and ( |28| ) have a 
solution and if it is unique. Of course only solutions such that Wq.^;) > are acceptable and these conditions introduce 
further restriction. 

In general, there are two possibilities. The first possibility is that there are no solutions which satisfy eqs. (p8|). In 
this case we can relax eq.([28|) by imposing 

E 4 i^^^h = E 4 1 iS^SJ)l I fc = 0, 1, 2..kM (30) 

ij ij 

where is the distance between spin i and j. Choosing in an appropriate way it is possible to reduce the number 
of conditions until a solution is found. If fc^f > 2 we believe that the solution is rather reasonable since the conditions 
of Sec. II are satisfied on the unit. 

The second possibility is that the solution is not unique. In this case we expect only small differences among 
different solutions. 

To show how this scheme works, in the following sections we will analyze a number of frustrated systems. In all 



cases we are always able to satisfy conditions (|16|) and (28) 



V. A DECORATED ISING MODEL. 

In order to check our approach we have considered a decorated Ising model with frustration. Starting with an Ising 
model on a square lattice we introduce between each pair of n.n. spins Si, Sj on the square lattice two extra spins 
S'fc, Si (see Fig. 5) modifying the interaction from 

n{s,,Sj) = -j{s,Sj-i) (31) 

to 

T-L{Si, Sj, Sk, Si) — —J{SiSj + SiSk + SiSi + SkSj ~ SiSj — 3) ; (32) 



this generic set of four spins will be the elementary unit I. For simplicity in eqs. ( |31| ) and ( |32| ) and in the following 
of this section we omit the label I. 

The partition function of this model is reducible to the Ising one via a "decimation" on spins Sk, Si: 
Z^Sfc Si e"^'^('^"'^J"^'=''^') = A(J)e^'^('^''^J~^^. The critical temperature can therefore be calculated exactly Tc = 2.24.... 
Monte Carlo estimation of the percolation temperature for the unmodified KF-CK clusters on such a system gives 
7p — Tc ~ 0.2 (here and in the following temperature are expressed in units of J/ks)- As expected the presence of 
frustration prevents the coincidence of percolation and thermodynamic properties for the unmodified KF-CK clusters. 



We have to solve eqs. (|16|) and (29) for the unknown weights Wa where a labels the bond configurations in the 
elementary unit (Fig. 5). The average in eq. ( ^ ) is over all spin and bond configurations with probability given by 
( p6| ) where Hi is given by (|3^). 



The spin correlations can be easily calculated from (27) since they does not require the knowledge of the Wa- We 
can immediately find (Fig. 5) 

{S.,Sj}i > (33) 
{S,Sk)i = {SkS,)i = -{S,Si)i = {SiS,)i > (34) 
{SkSi)i = 0. (35) 
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We can disregard bond configurations by inspection. For example the weight of a bond configuration which connects 
i and j through site / must be zero. In fact this bond configuration would correspond to Si and Sj antiparallel resulting 



in {'^fj)i > contrary to eq. (|29|). By imposing eqs. (29)-([35[) we reduce the number of possible bond configurations to 



12. Furthermore 3 of them have the same connectivity properties (i.e, they connect the same sites) as the configuration 
a — 3 (Fig. 6). Therefore they can be disregarded reducing the number of weights different from zero to 9. Eqs.(p^) 
now read 



Wq + W4 = 



Wq +Wl + W4 
Wo + U'2 + 3U'4 
Wo + Wi + W2 + W3 + 3W4 



(36) 



where u = e"^^''. The structure of such an equation system still allows us to chose one unknown arbitrary; in 
particular we have imposed — because it provides Wi > VT. 
Then the solution is 

Wo = u'^, Wl = u^{l - u), W2^U^{1-U^) 

W3 — u{l — U — + U^), W4 — . ^ ' 

In order to calculate the percolation temperature of the clusters defined by eq. (|37|) , Tp, we proceed in the following 
way. Given a spin configuration {Si} we assign to each plaquette I a bond configuration with the probability P{ca\{Si}) 
provided by eq. (^2|) and (^^ . Then we obtain clusters defined in the entire lattice as maximal sets of spins which are 
connected by bonds (CK-like cluster definition). It is, then, possible to measure percolation quantities. 

We have estimated Tp via a data collapsing of the probability P{T, L) of having a percolating cluster at temperature 
T in a system of size L (Fig. 7). We have simulated the model with both standard spin-flip MC dynamics and the 
cluster MC dynamics which we will discuss in sec. VIII. A obtaining indistinguishable results within our numerical 
precision. In Fig. 7 we have reported the results obtained with the latter dynamics. Using the scaling "ansatz" that 
near the percolation point P(T, L) = f{{T — Tp)L^/^p), where the functional shape of / is unknown, we have found 
that Tp ~ 2.25 and i^p ^ I consistent with the critical temperature Tc — 2.24 and with the critical exponent i/ — 1 
a ferromagnetic Ising model. 



VI. THE ASYMMETRIC FULLY FRUSTRATED MODEL (AFF). 



Let us consider a less trivial model, the frustrated Ising model on a square lattice with periodic boundary conditions 
where each plaquette contains three equal ferromagnetic interactions J and one antiferromagnetic interaction —XJ 
(0 < X < 1). This model interpolates between the Fully Frustrated (FF) model {X = 1) and the diluted Ferromagnetic 
Ising model (X = 0) |§. 

If we take Wa according to the definition given in Fig. 8 eqs.([l^) specified to this case give 

Wl + U>3 = 

wi + W2 = y?^-^ 

Wl + 2W2 + W3 + w^ + We + Wj + Wg — u (38) 
Wl + 2W2 + ^3 + U>4 + 2wj + Wg ~ U 

Wl + 3w2 + u>4 + 2we + ws = 

where u — e~^'^"^. The number of unknowns is larger than the number of equations. The percolation temperatures 
associated to these solutions will be, in general, higher than the thermodynamic critical temperature. In order to 
have these two temperature as close as possible and the cluster connectivity representing as better as possible the 
spin correlation, as discussed in sec. IV, we can impose either conditions ( psf ) or more simply eq. (p9|) . For pedagogical 
reasons we choose here eq. (ESh. These are three independent equations 



8 



Aw2 + 2w4 + 2wq + 2w7 + ws + Swg = - + u - (39) 
4w3 + 2w5 + 3w7 + ws + 3w7 = \u^ -u^ -3u + 3u'^+^\ 



Eqs.(g8D and (|39|) have the solution ^ 

Wi — U^, W2 = v?{u^ ~ u) (40) 

u>4 = u(l + - 2m-^+1) (41) 

[0 u < u j ^ ' 

-\2w, u<u*] 

_ / u> u* \ , , 

^8-| U^_u3_3^ + 3^X+2 y<y*| (44) 

and ^3 = uij = uig = 0. The equations ( p8| ) and ( |39| ) do not provide a unique solution. In fact a general solution can 
be found choosing as a free parameter. However, the further requirement Wa > leads to = 0. The solution 
changes form for a temperature T* such that u* = ^-"^-^I^bT satisfies the equation (1 + 3u^)u"^ — u'^ — 3u = 0. This 
is due to the fact that at T = T* the correlation between spins linked by the interaction ~XJ (let's call them 5*1 and 
^2) changes sign leading to different possible bond configurations. For example, the configuration associated to the 
weight W5 which links antiferromagnetically Si and S2 vanishes when {S1S2) > [^ . 

It is interesting to specify the general solution for the diluted ferromagnetic Ising model {X = 0) and the FF model 
(X = l). 

In the first case {X = 0) it results u < u* for any temperature and we get the following non zero weights 

wi = W2 — u^(l — u) W4 = We = u{l — u)^ ws = {1 — u)'^ . (45) 

These weights reproduce the original KF-CK solution for a plaquette; in fact they are in the form of products of u 
and 1 — u which are the KF-CK weights for a single spin pair. Since for X = the model reduces to a ferromagnetic 
Ising model where some interactions have been put equal to zero, the original KF-CK solutions reproduces the right 
clusters. 

In the case X = 1 it results u > u* for any temperature and one obtains the following non zero weights 

wi = u'^ W4 — ~ u{l — u^) (46) 

which are in agreement with the cluster structure used in Ref. |^ . It is worth while to note that in this limit all the 
bond configurations which connect spins on opposite corners have weights equal to zero preventing the four spins of 
the plaquette to belong to the same cluster even at T = 0. 

In order to check if the percolation model we have defined has the expected properties, i.e. if Tp = Tc, we have 
studied the percolation and spin properties of the system with both spin-flip standard MC dynamics and the cluster 
MC dynamics which we will discuss in sec. VIII. B obtaining indistinguishable results within our numerical precision. 
In Figs. 9 and 10 we have reported the results obtained with the latter MC dynamics averaging at least over 6 • 10''' 
MC sweeps after discarding the first lO''. 

In Fig. 9 we show the data collapsing for the mean cluster size defined as S' = rigs'^ where Hs is the number of 
clusters of size s and the sum is over finite clusters for three values of X {X = .5, X = .75, X ^ 1). For comparison 
in Fig. 10 we show, for the same X values, the data collapsing for x — ((-^^) ~ (l-^l)^) where M is the magnetization. 
From those data we extract the percolation temperature, Tp{X), the critical temperature Tc{X) and the critical 
exponents "fp(X), h'p(X), "/{X), and v{X) for any X value. Summarizing we find 
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Tp(0.5) ~ 1.24 rp(0.75) ~ 0.97 Tp(l) < 0.1 (47) 

7p(0.5) ~ 1.75 7p(0.75) ~ 1.75 7p(l) ~ 2.0 (48) 

iyp{0.5) ~ 1.0 i/p(0.75) ~ 1.0 i/p(l) ~ 1.0 (49) 

for the percolation quantities, and 

rc(0.5) ~ 1.24 TciO.75) ~ 0.972 T^il) ~ (50) 
7(0.5) ~ 1.75 7(0.75) ~ 1.75 7(1) ~ 1.51 (51) 
1/(0.5) ~ 1.0 j/(0.75) ~ 1.0 1/(1) ~ 1.0 (52) 

for the thermodynamic quantities. We note that Tp{X) 2± Tc{X) within the estimated numerical precision. Fur- 
thermore the same critical exponents control the divergence of percolation and spin properties 7p(^) — 7(X) and 
i^p{X) = i^{X) for all X values but X ^ 1. In the X = 1 case we find 7^ i^p{l) = 1.0 and 7(1) 7^ 7p(l) = 2.0. 
This result arises from the fact that the condition \ {SiSj) \ = {--fij) is satisfied only approximatively js^. The value of 
Tc{X), 7(X) and i'{X) are in agreement with the exact solution which gives rc(0.5) = 1.239..., rc(0.75) = 0.972..., 
Tc{l) = 0, Q J/ = 1 and 7 = 7/4 for X 7^ 1, ;/ = 1 and 7 = 3/2 for X ^ I. It is worth noting that different 
choices of the clusters which do not satisfy eq.(^^ as the unmodified FK-CK clusters systematically give Tp > and 



percolation critical exponents consistent with those obtained in the random percolation. |21 2 



A. Fractal structure. 



We have analyzed the cluster structure at T ~ Tp{X) ~ Tc(X) |20| , p5|| . For X ^ 1 we have found that a typical 
configuration of critical clusters is a fractal made of a backbone and dangling bonds. The backbone is made of links 



and blobs as found in the ferromagnetic Ising model |34| where the fractal dimension of the entire cluster was found 
to be equal to D = ^(^ + c^) which for dimension d = 2 gives D = 15/8 and the fractal dimension of the links or red 
bonds was found equal to Dr = 13/24. 

For the symmetric fully frustrated model X = 1 the structure changes drastically. In fact all the clusters are made 
of self-avoiding chains with fractal dimension given by the scaling relation D = ^(^ + 2). Using the numerical result 
7p ~ 2.0 and Vp ~ 2.0 we find numerically _D ~ 2 in agreement with the result of Ref.BADK,Coddington,Kerler 
that for T — Q predicts two percolating self-avoiding chains which fill up the entire system. It is interesting to note 
that if the approximation could be improved such that the condition (H) would be exactly satisfied we would obtain 
7 = 7p = 7/4 and v = Vp = 1. In the plausible event that the exact clusters are still self-avoiding chains the fractal 
dimension would be D = 7/4 identical to the fractal dimension of a self-avoiding random walk at the 9 point [5^. 



VII. ISING SPIN GLASS. 

The most complex and interesting model in the class of spin systems described by the Hamiltonian (|l|) corresponds 
to the case in which e^j is a random variable: this introduces quenched disorder together with frustration. We 
have studied the ±J Ising Spin Glass, in the case of a square lattice with the probability distribution p{tij) = 
kS f:-^^-i -|- (1 — k)S^.^^i where k is the concentration of antiferromagnetic interactions. The phase diagram of such 
system was described by Ozeki and exhibits at low temperature a paramagnetic-ferromagnetic transition if the 
concentration of antiferromagnetic interactions k is enough diluted (i.e. if < k < kq, with kq ~ 0.1), otherwise there 
is a spin-glass transition at zero temperature (i.e. ii kq < k < 0.5). The case k = 1/2 correspond to the Edwards and 
Anderson model (EA). Q 

As in the AFF case, we partition the lattice in square plaquettes of four spins. The system is then characterized by 
two kind of plaquettes: frustrated and unfrustrated. Frustrated plaquettes are those with one or three antiferromag- 
netic interactions |^^. We have analyzed such cases in section VI where we have obtained the probabilities P{ca\{Si}) 
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for the bond configuration, Cq, given a configuration of spin {Si} (see eq.(|2|) and (^)). We note that the weights 
given in eq. (|4^) for a frustrated plaquette containing three ferromagnetic and one antiferromagnetic interactions only 
depend on the interactions satisfied. Therefore the same weights can be used for the frustrated plaquette containing 
an odd number of antiferromagnetic interactions. We have used the unmodified KF-CK weights p ~ 1 — e^^'^'^ for 
the unfrustrated ones (eq.(p5|)). 

In order to compare the results obtained using the plaquette of four spins as a unit (we call this choice 2nd order 
approximation) with those obtained by using the single pair of spins as unit (1st order approximation) we have 
simulated the model with both standard spin-flip MC dynamics and the cluster MC dynamics which we will discuss in 
sec. VIII. C. The percolation quantities which we have measured, Tp, 7p and Vp, do not depend, within our numerical 
accuracy, on the dynamics used. 

In the case of EA model (k = 0.5), we have found a percolation temperature Tp ~ 1.40 higher than the critical 
temperature = 0, but lower than that obtained with the unmodified KF-CK clusters where Tp^^ ~ 1.80. [0 We 
have also estimated the percolation critical exponents Vp and 7p via a data collapsing (see Fig. 11), obtaining the values 
Vp ~ 1.33 and 7^ ~ 2.36, which are consistent with the random bond percolation values Vp — 4/3 and ^p = 43/18. 

We have also studied the region of low k where Tc is finite and the transition is ferromagnetic. In Fig. 12 we 
show Tp^\k), Tp'^\k) and Tc{k) for < k < 0.1 and n — 0.5. The values have been obtained after a data collapsing 
for the mean cluster size 5* (eq.(19)) and susceptibility Xj respectively. It's clear from Fig. 12 that the percolation 
temperature T^\n) is again lower then the one obtained for unmodified KF-CK clusters Tp^\k), and higher then the 
critical temperature Tc{n); for values k > 0.1 percolation temperatures slowly decrease reaching the k = 0.5 value: 
Tp'^^ (0.5) ~ 1.8 and T^^\o.h) ~ 1.4, while abruptly goes to zero. || 

From this analysis it comes out that neither the unmodified KF-CK clusters nor our clusters are able to correctly 
represent spin correlations in Spin Glass systems. However since the percolation temperature Tp^-* < Tp^^ one might 
expect that a systematic improvement can be obtained if larger elementary units are used (Fig. 2). 



VIII. MONTE CARLO DYNAMICS ASSOCIATED TO PERCOLATION MODELS 



It is possible to apply the general definition of cluster given above, to develop general MC cluster dynamics. The 
clusters are constructed assigning to each elementary unit one of the possible bond configurations according to the 
probability given in eq. (p2[) . Then the usual SW algorithm can be applied to the clusters described above. Following 
Ref. 01 it can also be proven that detailed balance holds (see Appendix). 



A. Decorated Ising modeL 



We partition the original square lattice in plaquettes as described in section V, and use the clusters defined there. 
With that cluster definition we have implemented the SW generalized cluster dynamics. 

We have estimated the percolation temperature Tp ~ 2.25 and the percolation critical exponents 7p ~ 1.77 and 
i/p ~ 1. The values obtained are indistinguishable from those obtained by using spin-flip MC dynamics reported in 
sec. V. We have also estimated the corresponding thermodynamic quantities Tc ~ 2.24, 7 ~ 1.78 and v ~ 1.05 which 
are consistent with a ferromagnetic Ising critical point within our numerical accuracy. 

In order to study the relaxation times we have computed the time dependent magnetization correlations 



{\ M{t)M{t')\) ~ {\M{t') \Y 
{M{t'Y) - {\M{t'W 



^ ,.,rf,i^2\ /I..^/.AI^2 (53) 



where M{t) is the magnetization at time t. Using the new cluster dynamic we find a dramatic reduction of the slowing 
down which is present for the standard SW dynamics. The new dynamics has critical autocorrelation times of about 
10 MCS (Monte Carlo Step per Spin), whose order of magnitude is comparable to those of the standard SW dynamics 
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for a ferromagnet of the same size at criticahty. On the contrary standard SW algorithm on Decorated Ising model 
shows very large correlation times at Tc (see Fig. 13). Our approach, then, reduces the critical slowing down in this 
system when compared to standard SW and local spin flip dynamics. 

B. AFF model. 

We partition the original lattice in elementary units as described in section VI and use bond configurations and the 
associated probabilities introduced there to define clusters. With such cluster definition we have then implemented 
a SW generalized cluster dynamics. We have estimated percolation quantities Tp(X), "tp{X) and Vp{X) for X = 
0.5,0.75, 1.0 obtaining values indistinguishable from those calculated by using spin flip MC dynamics. We have also 
computed the critical temperature Tc{X) and critical exponents ^{X) and i'{X). 

In order to study the relaxation times of our SW generalized cluster dynamics we have calculated the magnetization 
correlation function ( |53| ) versus Monte Carlo sweeps at the critical temperature Tc{X). It shows a dramatic reduction 
for all the X values which we have studied {X = 0.5,0.75,1.0) with respect to SW unmodified CK-KF cluster and 
local MC dynamics. For a quantitative analysis of the critical dynamic exponent z{X) at Tc{X), we have calculated 
the integrated autocorrelation time t using the self-consistent procedure suggested in Ref. | |39t with a window equal 
to 6. This method allows us to calculate r for different system size L in a consistent way. We found a power law 
scaling of the form t = kL^ as shown in Fig. 14. The estimated values of z{X) for X = 0.5,0.75 are z{0.b) ~ 0.30 
and z(0.75) ~ 0.46. The result definitely shows a strong systematic reduction of the critical dynamical exponent 
compared with those of standard MC dynamics. These results seem to indicate that the criterion to let Tp{X) as near 
as possible to Tc{X) (we have previously showed the coincidence of this two temperature for these models), allows to 
individuate efficient cluster dynamics. It is worth noting that even in the case X — \ when Tp = = Vp — i' — 1 
and 7p > 7 the cluster dynamics exhibits a drastic reduction of the critical slowing down. Nevertheless, our analysis 
suggests that, for X = 1, it is possible improve further this result considering larger units as starting point of the 
proposed procedure. 

C. ± J Ising SG model 

The panorama is more variegated in the more complex case of Ising Spin Glass with varying ferromagnetic inter- 
actions concentration k. This model exhibits in 2d a paramagnetic- ferromagnetic transition for Q < k < Q.l and a 
spin-glass transition for 0.1 < k < 0.5. Analogously to the other presented cases the cluster dynamic is realized by 
using the clusters defined in sec. VII. For frustrated plaquettes we have used the probabilities calculated in sec. VI for 
X = \ (eqs. (^6|) and (|2^)), while for unfrustrated plaquettes we have used unmodified KF-CK clusters (eqs. ( |45| ) 
and (H)). 

To check our simulation, we measured thermodynamic functions as energy E and specific heat Cy. pO[ | They 
reproduce known data in literature up to a temperature, Tf, under which our MC cluster dynamics freezes. We have 
also estimated Tp, 7^ and Vp obtaining good agreement with the values obtained by using spin flip MC dynamics. 

We have already noted that in Ising Spin Glasses the percolation temperature of unmodified KF-CK clusters, 
Tp^\k)^ is higher than the percolation temperature, T^\k), of the clusters defined in section VII: we have 

T«(«)>Tf (^)>r,(^) (54) 

where Tc is the thermodynamical critical temperature. The equality holds only for k ~ Q and 1 (ferromagnetic and 
antiferromagnetic case respectively). 

It is possible to summarize the results in this way: as n departs from the ferromagnetic Ising model k = 0, relaxation 
times for temperature close to the critical temperature Tc{k) get longer, and in the region where ferromagnetic phase 
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disappears (0.1 < k < 0.5), they become extremely long, even if always shorter than those of both standard SW 
cluster dynamics (unmodified KF-CK clusters) and local spin flip dynamics. 

Along the paramagnetic-ferromagnetic transition line (i.e. at Tc{k) with < k < 0.1) we have estimated the critical 
autocorrelation time t(k), defined as the time to reduce square magnetization correlation to 1/10 of its value at t — 0. 
These results are shown in Fig. 15 for a square lattice of size L = 32. 

In the region where the ferromagnetic phase disappears (0.1 < k < 0.5) and the SG transition takes over at Tsg = 0, 
our simulations get worse. We have studied for the case k = 0.5 the following relaxation function 

l(t) = ^J2(S.^(^o)S,^t + to)) ■ (55) 

i 

as a function of time (Monte Carlo Step) for systems whose size is L — 80, 90, 100. Due to very long autocorrelation 
times we were able to perform simulation up to T/ ~ 0.8. Averages in eq.(p5|) were taken over 1 ~ A ■ 10"* MCS 
discarding the first 5 • 10'^ MCS. We observed that relaxation time of our cluster dynamics above Tf is at least one 
order of magnitude lower than that of a standard spin flip dynamics (see Fig. 16 and, for comparison, Ref. |^). 

In conclusion in the case of spin glass we see that to a lowering of the difference \Tc — Tp\ corresponds a reduction 
of the relaxation times. However, there are indications that such a reduction exists only for T > Tf. Our results 
suggest that taking a partition of the lattice made by larger "elementary" units the procedure we have discussed 
define clusters whose percolation temperature is closer to the critical temperature of the original spin model. The 
associated cluster dynamics is expected to be characterized by shorter autocorrelation times. Work is in progress in 
this direction. 



IX. SUMMARY. 

In this paper we have discussed a general scheme for devising efficient MC cluster dynamics for spin models. The 
scheme is based on three main steps. The first one consists in choosing a partition of the lattice into "elementary 
units" . Then, using a method first introduced by Kandel, Ben-Av and Domany |l^ which is based on independent 
choices on each elementary unit, it is possible to define a vast class of cluster models whose Free Energy is identical 
to the original spin model. Finally, among the many cluster models it is possible to choose the one which satisfies 
at the best the equality between cluster connectivity and spin correlation (cfr. eqs.(^) and (|28|)). This procedure 
defines clusters which can be used to implement a MC cluster dynamics. We have applied this method to a number 
of 2d frustrated spin models taking as elementary unit a single plaquette. We show that every time Tc ~ Tp (i.e. the 
thermodynamic critical temperature of the spin model Tc is equal to the percolation temperature for the equivalent 
percolation model) the associated MC cluster dynamics is characterized by very small autocorrelation times and a 
critical dynamic exponent z much smaller than the one obtained in local (Metropolis) MC dynamics. We have also 
shown that in the more complex case of a spin glass where disorder is added to frustration the percolation temperature 
Tp results larger than spin glass temperature Tsg and the percolation critical exponent are consistent with random 
bond percolation exponents. However we see that in this case the percolation temperature can be decreased up to 
Tp ~ 1.4 using a lattice partition based on a single plaquette. This result suggests that taking larger elementary 
units, like the ones in Fig. 2d, Tp can be further reduced. Then, we still find a lowering of the autocorrelation time for 
T > Tf ~ 0.8 but, this time, there are no indication of a lower z compared to standard Metropolis dynamics. 

In conclusion our procedure allows for a systematic decrease of the autocorrelation times and, therefore, may serve 
as a general framework for the development of efficient MC dynamics in frustrated spin models. 



APPENDIX 



The aim of this appendix is to show that the MC dynamics defined in sec. VIII satisfies detailed balance. Following 
Rcf. jl^ we show that provided the mapping from Ti to ii (see eq. (|l6|) ) , a MC dynamic which verifies detailed balance 
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principle for Ti verifies it also for Ti. In particular after executing a freezing and deleting operation on the original 
spin system, we can implement with the built clusters a cluster dynamic based for example on random flipping of 
independent clusters as in SW procedure. This is possible because such a dynamics certainly satisfies the detailed 
balance principle and is generally ergodic at finite temperature. 

To prove that detailed balance is respected, let us make the following preliminary considerations. We can rewrite 
the relation (|6|) as 

a(l) 

This is the normalization condition for the conditioned probability to have the bond configuration c^^i) on the l-th 
elementary block, given the spin configuration {Si} on the system, i.e. the (^6|) expresses the normalization condition 
for the 



Wr,(l)e 



nc.(im}) = "'\_pn.{s.} • (57) 

Since the choices on the elementary block are independent, the probability of the bond configuration C on the whole 
system, given the spin configuration {Si}, is the product P{C\{Si}) = Hz P{'^a{i)\{Si})- 
To obtain detailed balance principle, we must impose the following condition: 

e-P^i{s^})T{{S,} -> {S[}) = e-'3«»^'»r({5a {S,}) (58) 

where T{{Si} {S'i}) is the transition probability from state {Si} to state {5*^'}. 
By definition T{{S'^} {Si}) may then be written as 

T{{S'i} ^ {S,}) = ^F(C|{S:})fc({S:} ^ {S.,}) (59) 
c 

where Tc{{S'^ {Si}) is the transition probability associated to the dynamic that we use on the dilute system with 
the Hamiltonian H. (for example this dynamic may be the simple random flipping of independent clusters). Generally 
let's suppose that the Tc({S'-} — > {Si}) respects the detailed balance principle, i.e. 



e 



'm{s[},c)f^(^{Sl} ^ {S.}) = e-''^«^'>^^)fc({5.} - {Si}) . (60) 



Therefore, from the (59) and the definition of P{C\{Si}), we have 



msl} - {S.}) = E ( n """^'^^ -0mis'} 1 ^ci{Sl} -> {S.}) , (61) 



-/3-H,({Sa,c„(i)) 



C 

that, using the (|6^), becomes 



T{{Sl} - {S,}) = E \-pn.iS'J j ^c:({5.} - {S',}) . (62) 



Now, multiplying and dividing for e we obtain 

mSl} {5.}) = -3^^^ [II ' ' ,-,n.is.} j Tc{{S.} - {SI}) , (63) 

and by ( |59| ) and the definition of P{C\{Si}), we get 

T{{Sl} {S,}) - _^^^^r({5,} ^ {5a) . (64) 
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This expression is the ( p8| ) and therefore the vahdity of the principle is demonstrated . 

Summarizing, the main assumptions underhng this proof consist in supposing that we can write the Hamiltonian 
as sum on elementary blocks Ti = 7i;, the choices on the elementary blocks are independent from each other and 
we are using a dynamics for the dilute system H, that respects the detailed balance principle. 

In particular a generalized cluster MC dynamics may be implemented with the following steps: individuate the 
clusters with "freezing and deleting" (i.e. to map Ti. into H); random flip of such clusters (this move certainly verifies 
detailed balance because clusters are not interacting in 7i), and iterate. 

About the ergodicity we can say that the cluster dynamics here described is certainly ergodic for every finite 
temperature, because the probability to go from a given spin configuration to any other is always different from 
zero for every non zero temperature. In general ergodicity at zero temperature is difficult to prove: it must be 
checked specifically in each particular case. Nevertheless it is possible to guarantee ergodicity also in such extreme 
conditions, alternating cluster moves with a dynamics that certainly is ergodic at this temperature, without changing 
the qualitative features of the cluster dynamics [Q . 

In conclusion we have proven that adopting the proposed mapping of Hamiltonians, and in particular for our general 
definition of clusters, it is possible to develop MC dynamics which verify detailed balance principle. 
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FIGURE CAPTIONS 

Fig.l A schematic example of two spin and bond configuration. The spin at site 1 belongs to the infinite cluster in 
both configurations (a,b) with different orientations. While both configurations give positive contributions to 
the percolation probability Pi they give opposite configurations to the magnetization (5*1). Similarly the spins 
at sites 1 and 4 are connected in both configurations (714 = 1). However, they are parallel (7I4 = 1) in a) and 
antiparallel {'yf^ = 1) in b). Therefore both configurations give a positive contribution to the pair connectness 
function pi4 but opposite contribution to the pair correlation function {S1S4). 

Fig. 2 Examples of possible elementary units partitioning a square lattice. The unmodified KF-CK clusters (see text) 
are constructed starting from elementary unit b). The clusters discussed in sec. V, VI and VII make use of 
elementary unit c). 

Fig. 3 A checkerboard partition of a square lattice. To cover all the lattice, we can choose the shaded plaquettes, or 
the unshaded ones, as set of elementary units. 

Fig. 4 The possible bond configurations on the square plaquette: full lines are infinite interactions or present bonds, 
while zero interactions are not marked. The label a is the index of configuration Ca and of the corresponding 
statistical weight Wa. As an example the configuration a = has no bond present. 

Fig. 5 The Decorated Ising model described in the text: each pair of interacting spin Si, Sj in a square lattice 
is decorated by spins Sk, Si. Full lines are ferromagnetic interactions, dashed lines are antiferromagnetic 
interactions. 

Fig. 6 The bond configurations of the decorated Ising model (sec. V) whose weights are different from zero (eq.57). 
We assign the same weight to all the elements belonging to the same group. In the figure each group is identified 
by a curl bracket. The conventional representation of interactions is the same as in the Fig. 5. 

Fig. 7 Decorated Ising model in two dimension (see sec. V): data collapsing for the probability P{T,L) of having a 
percolating cluster at temperature T in a systems of size L = 10, 20, 30 (for each size the number of system spins 
is 5 X L^). The data have been obtained by using the cluster dynamics discussed in sec. VIII. A. 

Fig. 8 Bond configurations for the AFF model (sec. VI). The 15 possible bond configurations are grouped, by symmetry, 
in 9 groups. To each element of a group is assigned the same weight Wa (cfr. eqs. (59) and (60)). In the figure 
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wc show only one clement for each group; the other elements can be obtained trivially conserving the number 
of ferromagnetic and antiferromagnetic bonds. In the brackets it is shown the number of bond configurations 
belonging to the specific group. Full (dashed) lines are bonds between ferromagnetically (antiferromagnetically) 
interacting spins. 

Fig. 9 The data collapsing for mean cluster size S of AFF models for systems with X = 0.5, 0.75, 1.0 and number of 
sizes L. The data have been obtained by using the cluster dynamics discussed in sec. VIII. B. 

Fig. 10 The data collapsing for susceptibility for the AFF model, with the same values of X and L as in Fig. 9. 

Fig. 11 The data collapsing for mean cluster size S for EA Spin Glass model for systems with size L = 48,64,80, 100. 
Critical exponents and temperature are also reported. 

Fig. 12 The critical temperature Tc(k) , the percolation temperature Tp^\k), for immodified FK-CK clusters, and the 
percolation temperature Tp {k) , for the clusters discussed in sec. VII, versus the antiferromagnetic interactions 
concentration k for an Ising Model with variable antiferromagnetic interaction concentration. The arrows show 
the values of fc = 0.5. The data reported have been extracted by data collapsing of mean cluster size, S, and 
susceptibility x (see text). The cluster dynamics used is described in sections VIII. C. 

Fig. 13 Correlation function, ^{t), as function of time (MC steps per spin) for the decorated Ising model introduced in 
sec. V at T = 1.6 and for a system size L = 32. Two MC dynamics are compared: unmodified FK-CK cluster 
dynamics (dashed line) and the cluster dynamics introduced in sections VIII. A (full line). 

Fig. 14 The relaxation time r versus system size L for the AFF model (sec. VI) with X = 0.5, 0.75. Assuming a power 
law scaling r = kL^ the estimated values of z{x) are z{0.5) = 0.30 and ^(0.75) = 0.46 

Fig. 15 The relaxation times r versus antiferromagnetic interactions concentration k at T = T^{k) for the Ising Model 
with variable antiferromagnetic interactions of size L = 32. T = T^{k) is defined as the temperature at which 
the susceptibility of a system of size L = 32 gets its maximum. Two MC dynamics are compared: unmodified 
FK-CK cluster dynamics (squares) and the cluster dynamics introduced in section VIII. C (triangles). 

Fig. 16 The relaxation function, q (see eq.79), for Ising SG model versus time (MC Steps per spin). The temperature 
reported are, from the bottom to the top, T = 1.4,1.3,1.2,1.1,1.0,0.9,0.8 and the system sizes are: L = 
80, 90, 100. The cluster dynamics used is described in sections VIII.C. 
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